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Spherically symmetric random walks in arbitrary dimension D can be described in 
terms of Gegenbauer (ultraspherical) polynomials. For example, Legendre polynomi- 



als can be used to represent the special case of two-dimensional spherically symmetric 



random walks. In general, there is a connection between orthogonal polynomials and 

b : 

semibounded one-dimensional random walks; such a random walk can be viewed as 
taking place on the set of integers n, n = 0, 1, 2, . . ., that index the polynomi- 
als. This connection allows one to express random-walk probabilities as weighted 
inner products of the polynomials. The correspondence between polynomials and 
random walks is exploited here to construct and analyze spherically symmetric ran- 
dom walks in D-dimensional space, where D is not restricted to be an integer. The 
weighted inner-product representation is used to calculate exact closed-form spatial 
and temporal moments of the probability distribution associated with the random 
walk. The polynomial representation of spherically symmetric random walks is also 
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used to calculate the two-point Green's function for a rotationally symmetric free 
scalar quantum field theory. 
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I. INTRODUCTION 



In two recent papers we proposed and analyzed a new kind of D- dimensional random 
walk that is well defined even when D is noninteger. This random walk takes place on 
a spherical lattice consisting of an infinite set of concentric nested spheres of radii P n , 
n = 1, 2, 3, .... We define region n to be the volume lying between P n _i and R n , with 
the central region, region 1 being the volume inside R\. If the random walker occupies 
region n at time t, then at time t + 1 the random walker must move out to region n + 1 
with probability P on t{ n ) or in to region n — 1 with probability Pi n (n). The probabilities of 
moving out and in are in proportion to the hyperspherical surface areas bounding region n. 
Let Sd{R) represent the surface area of a D-dimensional hypersphere 

2tt d / 2 



We then take (for n > 1) 



SD(Rn) _ Rn ' 

S D (R n ) + S D (R n ^) ~ RD-i + i^-i 



^ut(n) = c/D r, v 7: D ^ = t^h^t (i.i) 



and 



S D (R n _i) R 



D-1 
n-l 



p. („) - ^ v - n ~ i n o) 

m[) S D (R n )+S D (R n ^) R^ + R^- 1 J 

For the special case n = 1 we define 

P out (l) = 1, P in (l) = 0. (1.3) 

Note that probability is conserved because the total probability of the random walker moving 
out or in is unity: 

P ut(n)+P in (n) = l. (1.4) 

To describe a random walk on this lattice we introduce the notation C n ,t- m , which rep- 
resents the probability that a random walker, initially in region m at time t = 0, will be 
found in region n at time t. The probability C n ^ m satisfies the partial difference equation 



Cn,t;m — Piny** + l)Cn+l,t-l;m + -fout(w — l)Cn-l,i-l;m (w > 2), 



(1.5) 



Cl,t;m — -Pin(2)C2,t-l ;m , (1-6) 

and the initial condition 

The random walk described above has the advantage that the quantity C n ^- m is a mean- 
ingful probability for all real values of the spatial dimension D; that is, for all times t, the 
inequality 

< C*i,t,;m < 1 

holds. This result is in stark contrast with the random walk as it is conventionally defined on 
a hypercubic lattice . For example, on a D-dimensional hypercubic lattice, the probability 
that a random walker who is initially at the origin will again be found at the origin at 
t = 2 is 

c - 1 

which is greater than unity for D < ~; the probability that a random walker who is initially 
at the origin will again be found at the origin at t = 4 is 

6D-3 



8D 3 ' 

which is negative for D < |. 

In this paper we present an array of new results concerning D-dimensional random 
walks. Specifically, we consider random walks that are defined by the partial difference 
equations (|1.5|), ( |1.6| ) and (|1.7|). We show in Sec. || that there is a natural one-to-one 
correspondence between the probabilities C n ^- m that describe a random walk on a lattice 
consisting of regions n, n = 1, 2, 3, . . ., and a set of orthogonal polynomials {Qn-i(x)}, 
n = 1, 2, 3, .... This set of polynomials is uniquely determined by the functions Pm{n) and 



-Fout(^) in Eqs. (p..5|), (|1.6|) and (|1.7|) . There is a simple expression for C ntt - m in terms of these 
polynomials. In general, one can view a random walk on the regions nasa random sequence 
of raising and lowering operators applied to the set of polynomials {Q n -i(x)}. (Although not 
discussed in this paper, this correspondence between polynomials and random walks extends 
to multidimensional random walks and multi- index systems of orthonormal functions.) 

If we take evenly-spaced concentric spheres (R n — n), we find that for the special cases 
D = 0, 1, 2, the polynomials {Q n -i(x)} associated with P; n (n) and P ou t(n) in Eqs. ( |1.1D , 
( |1.2|) and ( |1.3|) are standard {||| classical polynomials: Gegenbauer polynomials for D = 0, 
Chebyshev polynomials for D — 1, and Legendre polynomials for D = 2. However, for 
all other values of D the polynomials have not been previously studied and are not found 
in any of the usual treatments of orthogonal polynomials. While we can generate these 
polynomials, we have not been able to determine their general mathematical properties, 
such as their weight function and interval of orthogonality. 

In Sec. [TTI| we modify the form of P m {n) and P ut( n ) in Eqs. (|1.1|), ( |1.2| ) and ( |1.3| ) by 



replacing these functions with their large-n asymptotic behaviors. The polynomials that 
we now obtain are well-known classical polynomials (ultraspherical polynomials) for all D. 
This allows us to find closed-form expressions for the probabilities C n j. m for all values of D. 

Taking the probabilities in Sec. [Ill] we then calculate in Sec. |V| extraordinarily simple, 
closed-form, analytic expressions for the probability of a random walker eventually returning 
to the region from which the walker started, the expected time for the walker to return 
to the initial region, and other space and time moments of the probability distribution 
C n .t ;m - (In contrast, in Ref. 0] after heavy analysis we were only able to obtain asymptotic 
approximations for these moments.) We also find that for integer D these moments exhibit 
the qualitative features (e.g. Polya's theorem) of random walks on D-dimensional hypercubic 
lattices. 

Our long-range objective in studying JJ-dimensional random walks is to understand 
critical behavior in quantum field theory. We would like to understand, for example, the 
transition that occurs when a self-interacting scalar 4 quantum field theory in space-time 



dimension D < 4 becomes a free quantum field theory for D > 4. One possible approach 
to such a problem would be to formulate a quantum field theory in terms of random walks 
|3||J . However, if we do so on a hypercubic lattice, it is not possible to study these random 
walks except for integer values of D, as we have discussed above. As a result, we cannot use 
a hypercubic lattice to examine the behavior of a quantum field theory near D = 4. Thus, 
we are motivated to investigate alternative kinds of random walks that may be consistently 
defined for all real D. Critical behavior has already been observed in a two-dimensional 
spherically symmetric random- walk model 0. In the next paper in this series || we study 
this critical behavior as a continuous function of D for all D > 0. We show in this paper 
that polymers adhering to D-dimensional curved surfaces exhibit critical behavior. 

Of course, a quantum field theory that is developed from a spherically symmetric random 
walk will itself be spherically symmetric. Such a theory is physically unacceptable because 
it violates causality. Nevertheless, the critical behavior that is observed in such a theory 
may well be a universal function of D and, at the very least, such a theory may provide 
some clues as to how a scalar quantum field theory can go from interacting to noninteracting 
at D = 4. In Sec. [V] we carry out some preliminary investigations of spherically symmetric 
quantum field theory. Specifically, we use the machinery of spherically symmetric random 
walks that is developed in Sees. [II [IV| to obtain the free two-point Green's function of a 
rotationally symmetric scalar quantum theory. 

II. CONNECTION BETWEEN POLYNOMIALS AND RANDOM WALKS 

In this section we propose and discuss the following quadrature solution to partial dif- 
ference equations of the type (|1.5|), Q1.6Q and (|IT7|): 

Cn,i;m = J ^ dx W {x)Q m _ 1 (x) , (2.1) 

where {Q n (x)}, n — 0, 1, 2, . . . , is a set of polynomials orthogonal with respect to w(x) 
on the interval — 1 < x < 1 and {v n }, n = 0, 1, 2, ... ,isa sequence of positive numbers. 
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The form of (|2.1|) incorporates the initial condition (|1.7|) in a natural way. We simply 
choose to normalize the set of polynomials {Q n (x)} so that 

/ dxw{x)Q n {x)Q m {x) = —8 n>m . (2.2) 
J -i v n 

With this choice of normalization we see that at t = ( |1.7|) follows immediately from the 
above statement of orthogonality: 

C n fl- m = v n -i J i dxw(x)Q n -i(x)Q m ^i(x) = S n>m . 

We now demand that the set of polynomials {Q n (x)} obey the recursion relation 

Pin{n + l)v n Q n {x) = v n „ 1 xQ n _ l (x) - P mt (n-l)v n „ 2 Qn-2{x) (n>2), (2.3) 
and the initial conditions 

Qo(x) = 1 (2.4) 

and 



The partial difference equations ( pTo] ), (|1.6| ) and ( |TT7|) for the probabilities C n ^ m are auto- 
matically satisfied so long as Eqs. Q2.3|) , ( |2.4|) and ( p.5|) hold. 
We will now show that 

Qn(l) = 1 (2-6) 

for all n > 0. This interesting property is a consequence of the conservation of probability; 
namely, the probability of finding the random walker somewhere on the lattice at an arbitrary 
time t is unity: 

oo 
n=l 

To establish Eq. ( |2.7| ) we merely sum Eqs. (p..5|), ( |1.6| ) and ( |1.7| ) over all n > 1, using Eqs. 
(IOD, (|1.4|), and (|1.7|). Assuming that the sum 



f( X ) = H V nQni 



X) 



?1=0 

exists in the space of distributions, we substitute the expression for C n 4;m in Eq. (|2.1| ) into 
Eq. ( gj7|) to obtain 

J dxty(x)a;*Q Tri _ 1 (x)/(a;) = 1. (2.8) 

Next, we compute the function f(x) directly from the recursion relation (|2.3|) by summing 
over all n > 2, using Eqs. ( |2.4j ) and ( |2.5|) . We obtain the following equation for f(x): 



(1 - x)f(x) = 0. 
The solution to this equation is a generalized function: 

f(x) = a5(x-l), (2.9) 



where 5(s) is the Dirac delta function and a is a constant. Substituting Eq. ( |2.9| ) into Eq. 
( p.8| ) gives the condition that Q m _i(l) is a constant independent of m for all m > 1. Finally, 
from the conditions ( |2.4j ) and (|2.5|) , we conclude that this constant is 1, 

Q m (l) = 1, (2.10) 

and we therefore obtain Eq. ( |2.6| ). 

The result ( |2.6|) enables us to find a simple formula for the set of numbers {v n }. We let 



x — 1 in the recursion relation ( |2.3|) to obtain 

P ia (n + l)u n = u n _i - P ut(n - l)v n - 2 (n>2) (2.11) 
and in the initial condition (|2.5| ) to obtain 

vi = — "T - :• (2.12) 



The unique solution to Eq. ( j2.11| ) that satisfies Eq. ( |2.12| ) is 



V ^ V °R^TT) (n> - 1] - (2 - 13) 



The value of v o is determined from the orthogonality condition ( |2.2|) at n = m = and 
the initial condition ( |2.4j ): 

1 



J_ x dxw(x) 

The result in Eq. ( 2.13| ) can be used to eliminate the numbers v n from the recursion 



relation Q2.3|) giving a much simpler recursion relation for the polynomials Q n (x): 

Pout(n)Q n {x) = xQn-^x) - P in (n)g„_ 2 (x) (n > 2). (2.14) 
The initial conditions in Eqs. (|2.4j) and ( |2.5|) also become much simpler: 

Qo(X) = 1, = x. 

This recursion relation generates polynomials that exhibit parity symmetry; that is, 
even-index polynomials are even functions and odd-index polynomials are odd functions: 
Q n {— x) = (—l) n Q n (x). From the orthogonality condition (|2.2|) one can then deduce that 
the weight function w(x) is an even function of x. As a consequence, we see from the integral 
representation in Eq. fl2.1|) that an n versus t table of values of the probabilities C„ )t;m has 
a checkerboard pattern with nonzero entries alternating with zero entries in both the n and 
t directions. Evidently, a random walker starting from the site m at t = can only reach a 
site n at time t if n + m + t is even. This parity condition is a consequence of the original 
definition of our random walk in which the walker must move in or out on every step and 
may not remain in the same region. 

It is interesting to examine some special cases of the polynomial solution for C n ^ m m 
Eq. ( |2.1D . We consider the case of equally spaced spherical shells R n = n and look at some 
particular values of the dimension D with P out (n) and Pin(n) given in Eqs. ([Lip, ( |1.2| ) and 

(ED- 

A. Special case R n = n, D = 1 

Here, 
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PoutW = ~ and P in (n) = ^ (n > 2) (2.15) 

and 

P out (l) = l and P in (l) = 0. (2.16) 
The polynomials {Q n {x)} are the standard Chebyshev polynomials of the first kind H: 

7q (a;) = 1, 7l(x) = x, ^(x) = 2x 2 — 1, 

T 3 (x) = Ax 3 - 3x, TAx) = 8x 4 - 8x 2 + 1, 



and so on. For these polynomials, iy(ar) = — x 2 and u„ = 2/7r, n > 1, and t>o = l/V. 

The random walk probabilities in Eq. (|2.1|) are given by 



2 r 1 1 

C„ )tim = - / dx =a: f 7;_i(x)7^_i(a:) (n > 2), 

vr j-i vl — x 2 



If 1 1 

Ci,t;m = - / dx =a: f 7^_i(x), 

7T 7-1 Vl — X 2 

which for m = 1 reduces to the particular solution 

^n,n+2,-l;l ~ jj^ + j _ 1 ) !2 n+2i-2 ^ " ^' 



r (2*)! 



t!t!2 2 *' 



given in Ref. ffl. 



B. Special case 7? n = n, D = 2 

Here, 

Ti 7?, — 1 

Pout(n) = and P in (n) = . (2.17) 

v J 2n-l K ' 2n- 1 V 7 

The polynomials {Q n {x)} are the standard Legendre polynomials f|: 
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V (x) = l, V 1 (x)=x, V 2 (x) = i(3x 2 - 1), 

V z (x) = -{5x 3 - 3x), V±{x) = i(35x 4 - 30x 2 + 3), 
2 8 

and so on. For these polynomials, w(x) = 1 and v n = 2n + 1. Thus, the random walk 
probabilities in Eq. are given by 



C„ jt;m = (2n - 1) y dxx^n-^Vm-iix), 



which for m — 1 reduces to the particular solution 

(2n-l)(n + 2j-l)! 
"•"• 2/ ;:l j!(2n + 2j-l)!!# ' 

given in Ref. Jl|. 



C. Special case R n = n, D = 



Here, 



TL — 1 Ti 

PoutM = and Pin(n) = (n > 2) 

v ; 2n-l 2n-l V ~ ; 



and 



P out (l) = 1 and P in (l) = 0. 
The first few polynomials in the set {Q n {x)} are 

Q (x) = l, Qi{x)—x, Q 2 (x) = 3x 2 - 2, 

QJ X ) = hl5x 3 - 13a;), Q 4 (» = i(105a; 4 - 115a; 2 + 16), (2.18) 
2 6 

and so on. For these polynomials we have chosen v n = 3(2n + l)/[4n(n + 1)], n > 1, and 
t>o = 3/4. The random walk probabilities in Eq. (|2.1|) are then given by 



C n ,t ; m = T-r — TT / dx w(x)x t Q n ^ 1 (x)Q m ^ l (x) (n > 2), 
4n(n — 1) J-i 
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3 f L 

Ci,t ; m = 4 _/ rfx«;(a;)z*(5 m _i(a;). 

The polynomials Eq. (|2.18|) are closely related to the standard Gegenbauer (ultraspheri- 
cal) polynomials, {C^)(x)}, with upper index a = 3/2 ||. These particular Gegenbauer 
polynomials satisfy the recursion relation 

in + 1)C%1? (x) = (2n + 3)a< 3/2) (x) - (n + 2)c£_ 7 ? (z) (n > 0) 

and the initial conditions Cq (x) = 1 and C± (x) = 3x. These Gegenbauer polynomials 
are orthogonal on the interval — 1 < x < 1 with respect to the weight function w(x) = 
1 — x 2 . The polynomial Q n+ i(x) satisfies the same recursion relation as these Gegenbauer 
polynomial C^^ 2 \x). However, it is generated from different initial conditions. We have been 
able to show that the weight function w(x) with respect to which the set of polynomials 
{Q n (x)} is orthogonal satisfies the integral equation 

i-i 1-xt 2 (1 -x)[ln(l + y/x) -ln(l - y/x)]' 
We do not know a closed-form solution to this equation. 



D. Special case R n = n, D = 3 

Now, 

n 2 (n — 1 V 

PnJn) = — and RJn) - 



2n 2 -2n + l v ' 2n 2 -2n + l' 

For this case we can calculate any finite number of polynomials {Q n {x)}'- 

Qo(x) = 1, Qi(x) = x, Q 2 (x) = ^(5x 2 - 1), 

Qs(x) = 7^(65x 3 - 29a;), Q A {x) = ^(1625x 4 - 1130a; 2 + 81), 

and so on. These polynomials are orthogonal and they satisfy the normalization constraint 
( p.lOj ). However, for this value of D (and for all values of D other than D — 0, 1, 2) 
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these polynomials are not related to the standard classical polynomials that one can find 
in reference books. We are unable to determine analytically the weight function w(x) with 
respect to which these polynomials are orthogonal! Thus, the formal expression Eq. ( [2 . 1| ) 
for the probabilities C n ^ m is not very useful. 

In the next section we devise a random walk process for which we can determine the 
weight function and thus find in closed form physically realistic probabilities C n ^ m for all 
values of D > 0. 



III. RANDOM WALKS FOR ULTRASPHERICAL POLYNOMIALS 

In this section we show how to modify the expressions for P ou t{n) and Pi n {n) in Eqs. ( |1.1D , 
( |1 . 2|) and ( |1.3|) so that we are able to obtain analytic closed-form expressions for C n ^ m for 
all values of D > for the case of evenly-spaced spherical shells R n = n. The random walk 
process examined in Subsection [II D| is too difficult to solve in closed form simply because the 
formulas for P ou t(^) and Pi n (n) in Eqs. (|1.1|) , ( |IT2|) and ( |1.3j ) become much too complicated 
when D takes on values other than 0, 1, or 2. 

As we will see, the polynomials generated by the recursion relation ( |2.14| ) belong to a 



set of well-known classical polynomials if we take the formulas for P ou t{n) and P m (n) to be 
bilinear functions of n of the general form 

p i \ an + b fo u 

Pout W = — , 3.1 

cn + a 

with Pi n (n) = 1 — Pout(^)- Note that bilinear functions contain three arbitrary parameters. 
We fix these parameters as follows. First, we demand that the random walk be confined to 
the values of n > 1. To impose this condition we require that P out (l) = 1, or, equivalently 
that Pm(l) = 0. This fixes one parameter. Second, we demand that the large-n asymptotic 
behavior of P ou t(^) in Eqs. ( |1.1|) and (|3.1|) agree to order n. These two conditions above 
yield the unique choice 

, ^ n + D — 2 , , . n — 1 , 

Pout in) = and P in (n) = . 3.2 

ouU ; 2n + P>-3 mK J 2n + D-3 y 1 
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By determining the arbitrary parameters in Eq. ( |3.1| ) at the two boundary points n = 1 
and n = oo we obtain a uniformly accurate approximation to P ut{ n ) and Pm{n) in Eqs. 
Q, O and flTD for all n > 1. In fact, Eq. Q agrees exactly with Eqs. Q), Q) 



and Q) for L> = 1 [see Eqs. ( pT5l) and (gTg )] and D = 2 [see Eq. ([TT7])]. For other values 



of .D, Eq. ( |3.2| ) continues to be a good approximation, as is verified in Fig. |I|, where we 
compare Pm{n) in Eq. (|1.2| ) with Pi n (n) in Eq. (|3.2|) for several values of D. 

The requirement that Eqs. ( |1 . 1|) and ( |3.1|) agree to order n as n — > oo incorporates the 
crucial dependence upon the dimension D of space; to wit, as D increases, a random walker 



is more likely to move outward than to move inward. As we will see in Sec. [TV], it is this 
bias that gives rise to Polya's theorem; this theorem states that for < D < 2 a random 
walker returns to the starting point with probability 1, while for D > 2 this probability is 
less than 1. 

Substituting the formulas above into Eq. (|2.14|) gives the recursion relation 



(n + D-2)Q n (x) = (2n + D - 3)xQ n _ 1 (x) - (n - l)Q n _ 2 {x) (n > 2). (3.3) 
Taking as initial conditions 

QoO) = 1, 
Qi(x) = x, 

we can easily use Eq. (|3.3|) to generate subsequent polynomials: 

g 2 (x) = i[( J D + i)x 2 -i], 

Q 3 (x) = i[(D + 3)x 3 -3x], 
Qi{x) = D 2 l +2 D ^ D2 + 8D + 15)x4 ~ (6jD + 18)x2 + 3 1' 

Q 5 (x) = D 2 1 +2 d ^ D2 + 12D + 35 ^ 5 ~ ( 10D + 5 °^ 3 + 15 ^' 
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These polynomials are just the Gegenbauer (ultraspherical) polynomials || Cf^'(x) normal- 
ized so that Q n (l) = 1: 



g.(x) = J! lr ^- 1 U ^. 



(3.4) 



r(n + £>- 1) 

Gegenbauer polynomials are hypergeometric functions in the variable (x— l)/2; furthermore, 
since they are polynomials, the Taylor series for Q n (x) about x = 1 terminates: 
q ( x ) = V M r(j + L> + n - l)r(D/2) fx - W 



^Viy rp + n-i)r(j + i)/2) V 2 '■ (3 ' 5) 

From the conventional theory of Gegenbauer polynomials || we immediately know the 
weight function with respect to which the polynomials Q n (x) are orthogonal: 



w 



(x) = (l-x 2 Y D - 2 ^ 2 . 



(3.6) 



Also, the normalization coefficients v n in Eq. (|2.2|) are identified as 

(2n + D - l)r(n + -D - 1)T[(D + l)/2] 



(3.7) 



v^Fn!r(£)/2)r(£)) 

Finally, we note that the polynomials Q n (x) satisfy the Sturm-Liouville eigenvalue differen- 
tial equation 



d 2 d 
'1 - x 2 )— - Dx^- + n(n + D-V 



doc doc 
and the first-order difference-differential equation 



Qn{x) = 



1 — x J - — V nx 

dx 



Q n {x) = nQ n . 



(3.8) 



Now that we have identified explicitly the polynomials Q n (x), the weight function w(x), 
and the normalization coefficients v n , we can use the formula in Eq. (|2.1|) to calculate the 
moments of C n ^ m and obtain a physical description of our random walk. 



IV. QUANTITATIVE DESCRIPTION OF THE RANDOM WALK 



In this section we discuss the properties of the hyperspherical random walk introduced 
in Sec. |T|. We calculate the probability of eventually returning to the starting point of a 
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random walk, the expected time of return, and various other moments of the random walk 
probabilities. As will be evident, the key advantage of this random walk is that all of these 
quantities can be calculated in closed form. 

A. Probability of eventual return 

In a physical description of a random walk the simplest and most natural question to 
ask is, what is the probability of eventually returning to the starting point? The probability 
that a random walker will eventually return to region m, given that the walker started 
in region m, is denoted U m (D). To calculate H m (D) we use generating function methods 
previously described [see Eq. (2.11) of Ref. 0]; to wit, 

Km(D) = 1 - — 5-i . (4.1) 

Our problem is now to evaluate the sum, which we denote by S m , in Eq. ( |4.1| ) . Using 
Eq. ( |2.1|) with Q n (x) given in Eq. Q3.4Q we have 

OO CO -i 

S m = J2 C m,2t; m = Y,Vm-l / dx (1 - [Q m _ X (x)] 2 . (4.2) 

t=0 t=0 J 1 

Note that this sum is divergent unless D > 2. To verify this assertion we observe that 
the large-t asymptotic behavior of the integral in Eq. ( |4.2|) is 

J 1 dx (1 - x 2 i D - 2 V 2 x 2t [Q m -i(x)} 2 ~ T(D/2)r D ^ 2 (t -> oo). 

Thus, 

U m (D) = l (0<D<2). (4.3) 

When D > 2 the sum S m converges, and we begin the evaluation by interchanging the 
order of summation and integration: 

S m = iw-i J\ dx (1 - x 2 )^/ 2 \Q m -i{x)} 2 . 
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We evaluate this integral exactly using the recursion relation Q3.3|) and the difference- 
differential equation ( |3.8| ). The result is 

Substituting into Eq. ( (4.1[ ) gives 



This result is exact for all m and Z) @. The probability in Eqs. ( |4.3|) and ( |4.5|) confirms that 
Polya's theorem holds for this model of a random walk, regardless of the region in which 
the walk begins. 

For a hypercubic lattice the probability of eventually returning to the starting point of 
a random walk is given in terms of an integral [see Eq. (2.12) of Ref. 0]: 

H(jD) = 1 ~ f~dte-<[I (t/D)]°' 

where Iq(x) is the modified Bessel function. Unlike the random walk discussed in this paper, 
when D > 2, n(D) cannot be given in closed form (except for the special case D = 3). 
However, the asymptotic expansion of 11(D) for large D is known [see Eq. (2.15) of Ref. [p]]]: 

U{D) ~ W + Wi + ■ ■ ■ <^°°>- 

Note that for large D, the probability of returning to the starting point of a random walk falls 
off algebraically like 1/D in both models. In contrast, for the hyperspherical random walk 
discussed in Ref. |l], the probability function Ui(D) = 1 — 1/((D — 1) falls off exponentially 
like 2 1 ~ D for large D. As functions of D, the hypercubic 11(D) and the hyperspherical 
H m (D) discussed here both exhibit cusps at D — 2. 

Observe that for large m, U m (D) approaches 1. This happens because the available 
entropy for the random walk becomes constant; at large radius a sphere looks locally like a 
plane. Indeed, as n — >• oo, the recursion relation (|3.3|) approaches that of a one-dimensional 



random walk for which P ut( n ) — Pm( n ) — h 



2 ' 
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B. Expected time of return 



As explained in Ref. [0], the expected time of return T m (D) of a random walker who 
begins the walk in region m is obtained from the first moment of C m ^ m : 



T m (D) = t=0 ^ m ' 2t ' m . (4.6) 




Again, using formulas (|3.3| ) and (|3.8| ) we can calculate the sums in Eq. (|4.6| ) straightfor- 
wardly. We find that 



Note that as D increases, T m (D) approaches 2, independent of the starting point m. This 
is because for very large dimension D, if a random walker does not return to the starting 
point on the second step, the random walker will never return; as D — > oo the entropy for 
moving outward dominates the walk. However, for fixed D as m increases T m (D) diverges. 
This is because for large m the D-dimensional walk approaches a one-dimensional walk for 
which the expected time of return is infinite. 



In general, all temporal moments can be calculated in closed form. The pth temporal 
moment J2uLo ^ p C m ,2t;m is a rational function of D and m whose complexity increases with 
p. The sum defining the pth temporal moment converges when D > 2p + 2 and diverges 
when D < 2p + 2. We list the first four temporal moments below. [Note that the zeroth 
moment, S m , is already given in Eq. (^.4|).] 



(0 < D < 4), 



T m (D) 



2(D-2)[(2m-l)D+2(m-l)(m-2)] 
(2m-l)D(D-4) 



(£>>4). 



C. Higher Temporal Moments 



t=o 



C m ,2t;m — S, 



oo 



(2m + D — 3)M ( 
D-2 




oo 



(2m + D — 3)M 1 
(D-A){D - 2)£>' 
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where 



~ 2 (2m + D-3)M 2 

2^ ^m,2t;m (£, _ 6 ) _ 4) _ 2 ) J D( J D + 2) ' 



- 3 (2m + D-3)M 3 

t=0 



(D - 8)(D - 6)(D - A)(D - 2)D(D + 2){D + 4)' 



M = 1, 

Mi = (2m - 1)£> + 2(m - l)(m - 2), 

M 2 = (2m - 1)L> 3 + 2(7m 2 - 13m + 7)D 2 + 4(m - l)(6m 2 - 17m + 16)D 
+12(m - l)(m - 2)(m 2 - 3m + 4), 

M 3 = (2m - l)D 5 + 2(19m 2 - 29m + 15)L> 4 + 2(96m 3 - 296m 2 + 386m - 173)£> 3 
+4(99m 4 - 486m 3 + 1093m 2 - 1184m + 477)£> 2 
+4(90m 5 - 621m 4 + 2040m 3 - 3683m 2 + 3414m - 1252)£> 
+24(m - l)(m - 2) (5m 4 - 30m 3 + 97m 2 - 156m + 104). 



D. Spatial Moments 

The kth spatial moment of a random walk is defined as a weighted average over the 
probabilities C nAm - 

oo 

(R% = J2n k C n , t . !m . (4.7) 

n=l 

Note that in general {R k )t depends on the starting point m of the random walk. We have 
suppressed the argument m because, as we will see, the leading asymptotic behavior of 
(R k )t as t — > oo and the first correction to this behavior is independent of m. (The second 
correction does depend on m.) 

We have found an exact expression for {R k )t for all values of t for the special case m = 1: 

(R k ) 2t = (2t + l) k + ]T X (-l) r+ 7(2r - 1, 2s - 1, 2t, k) (4.8) 

r=l s=r+l 
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and 



t t+i 



where 



(R k ) 2t+1 = (2t + 2) k + £ E ("l) r+ 7(2r, 2a, 2t + 1, fc), (4.9) 

r=l s=r+l 



ffj ,. f M- (x-y)r(o; + y -l)r(t+l)(x fc -^) 

2*^(x)^( ^ /)^(^)^(^)^(^±2)( J D + x + Z /-3)■ l ' ; 



This formula has the virtue that the D dependence is very simple; the parameter D 
occurs just once in the denominator of / in Eq. ( |4.10| ). Furthermore, for the special case 
of the zeroth moment, setting k = in Eqs. ( |4.8| ) or ( |4.9| ) immediately gives the result 
(R°)t = 1, which states that probability is conserved. For k > 1 this formula is inherently 
complicated. It is not easy to determine the asymptotic behavior of (R k )t for large t from 
Eqs. (|4.S| ) or ( |4.9|) because terms in the double sum oscillate in sign. 

To find the asymptotic behavior of (R k )t as t — > oo we use generating-function techniques. 
We rewrite Eq. ( |4.7|) as a derivative operator applied k times to a power series: 



J OO 

{R k ) t = \im(z—) k £ z n C n , t , m . (4.11) 

dz n=l 

Next, we substitute into Eq. ( f4.11| ) the integral representation for the probability C njt;TO in 



Eq. Q) and use Eq. We obtain 



= lim f 1 dx (1 - x 2 )(D - 2)/2x t Q m _ 1 (x)( 



oo 

z z n v n Q n (x) 

. n=0 



(4.12) 



It is convenient to use the expression for v n in Eq. (|3.7| ) and the recursion relation (|3.3|) for 
Q n (x) to evaluate the sum in Eq. (|4.12|) : 



E ^VnQnix) = ,\ \L (l - Z 2 )(l - 2xZ + Z 



2\-(D+l)/2 



We are interested in the behavior of the resulting integral as t — > oo. By Laplace's 
method this integral is dominated by values of x near 1 in this limit. Thus, for fixed m 
Eq. Q2.10D implies that we may replace Q m -\{x) by 1 to leading order; we thus conclude 
that the leading asymptotic behavior of (R k )t is independent of m. To obtain higher-order 
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terms in the asymptotic expansion we replace Q m -i(x) by the expansion in Eq. (|3.5|) . A 
straightforward asymptotic analysis of the resulting integral gives the first few terms in the 
asymptotic expansion of (R k )t for large t with m fixed: 

' D+fc-T 



(m - \){m + D - 2) 



D 

{21)- 1 + O (V 3/2 ) ) (t->oo). (4.13) 



2 / v V 2 

3L> 2 £; - 18DA; + 12D - 2k 2 + 33k - 28 
+ 12(D + k - 2) 

Observe that the leading term in this asymptotic expansion is precisely the same as the 
result in Eq. (3.4) of Ref. [§ for the case of spherically symmetric random walks described 
by the probabilities in Eqs. (|1.1| ), Ql.2|) and (|1.3| ) with R n = n. The result in Eq. Q4.13|) 



is obtained directly and with considerably less effort than that in Ref. 0, where only the 
leading asymptotic behavior was obtained. Note that the first two terms in Eq. |4.13| are 
independent of the starting point m. To verify the accuracy of this asymptotic expansion 
we compare the first three partial sums of this series with the exact values of the moments 
obtained numerically at t — 1000; this comparison is given in Tables | and ||. In Table | we 
consider the fcth moment for various values of k and D with m = 1. In Table |I| we consider 
the first and second moments for various values of the starting point m with D = 2. 

From the asymptotic behavior in Eq. ( |4.13 ) with k = 2 we can determine the Hausdorff 



dimension Dh of the random walk 0. We find that 

D H = 2 

for all values of D. This result agrees with that obtained in for a .D- dimensional hyper- 
cubic lattice. 



V. APPLICATION TO QUANTUM FIELD THEORY 

One of our long-range goals in our study of D-dimensional random walks is a deeper 
understanding of D- dimensional quantum field theory. In particular, we are interested in 
how critical phenomena in such theories depend on the dimension of space-time. We are 
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especially interested in how a 4 scalar field theory becomes free as D — > 4. We have 
already conducted several investigations of D-dimensional quantum-mechanical and field- 



theoretic systems fllO|-|T3|1. In this section, as an elementary illustration of how to apply 
our work on D-dimensional random walks to quantum field theory, we use the random 
walk probabilities C n ^- m in Eq. ( |2.1| ) with the polynomials Q n {x) given in Eq. ( |3.4|) to 
calculate the Euclidean two-point Green's function of a D-dimensional free scalar quantum 
field theory having spherical symmetry. We will then verify our calculation by taking the 
spherical average of the two-point Green's function of a conventional translationally invariant 
(nonspherically symmetric) Euclidean field theory [|Hjj . 



A. Derivation of spherically symmetric propagator from random walk probabilities 

Cn,t;m 

For this calculation we follow the standard recipe discussed in Ref. ||. Specifically, we 
begin with the generating function G(n, m, A) for the temporal moments of the probabilities 

C n ,t;m- 

oo 

G(n,m, A) = J^^n^m- 

t=0 

Our objective is to find the continuum limit of this expression. 

For definiteness we choose n, m, and t to be even: n = 2N, m = 2M, t = IT . Also, 
without loss of generality, we take N > M. Substituting the formula for C m ^ m in Eq. (|2.1|) 
with Q n (x) given by ultraspherical polynomials in Eq. ( |3.4j ), and w(x) in Eq. ( |3.6| ), and v n 
in Eq. ( |3.7| ), we obtain 



G(2A,2M,A) = -^±± dx{l-x 2 )^- £ (x\) 2T C { 2N 2 J (x)C { 2M 2 J (x). (5.1) 

nM u Jo t=n-m 

Next, we perform the sum in Eq. ( |5.1| ): 

G(2A,2M,A) = -jA-l I jf dx (1 - x 2 )^ ^ x2x2 C&l(x)(£j:}(x). (5.2) 
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To prepare for taking the continuum limit we make use of the equivalence of Gegenbauer 
and Jacobi polynomials [|J: 

^ _ r(D/2)T(2N-2 + D) {2? t B? )( . ( . 

W - _ _ 1 + D/2) > 2N-1 W- [5.6) 

Substituting Eq. ( ^.3| ) into Eq. ( |5T2"| ) and taking A" and M large gives 

(xX) 2N ~ 2M 



G(2N, 2M, A) = 2 i - D N D / 2 M l ~ D ' 2 cfe (1 - x 2 )^— 

Jo 1 



x 2 A 2 



x P^ii 2 } (x)^ M 2 -i 2 \x). (5-4) 

When A" and M are large the integral in Eq. (|5.4j) is dominated by values of x near 1. 
Thus, we make the change of variable x — 1 — e 2 s 2 /2, where e is a small parameter: 

G(2iV, 2M, A) = 4 £ 2 JV^'-^ £** * ''^v^ ^"2" 

x (f )^-> (1 - (f )^K-> (. - ■ <-) 

We now make use of the following asymptotic limit for Jacobi polynomials [§J: 



where J a {s) is a Bessel function. Because iV and M are large, we can use this asymptotic 
limit twice in Eq. ( |5.5| ): 

rsPilt — c 2 c 2 /n\-\2N-2M 

G(2N, 2M, A) = Ae 2 N D ' 2 M^ D l 2 / ds ( 7 A ^2 

Jo 1 — A 2 (l — e A s A /2y 

x jD_2[e(2N -l)s]jD_2[e(2M - l)s). (5.6) 

We introduce the continuum variables r and r' by 

/xr = e(2iV-l), ^ = e(2M-l), 

where // is a mass parameter. Note that r > r'. Also, since e << 1, we may replace the 
upper limit of integration in Eq. ( |5.6|) by oo and simplify the integrand: 

G(2N, 2M, A) = 2er D / 2 (r') 1 - D/2 /i / ^ r^- ^J^rs) Jn^(fxr's). 

Jo 1 — A / (l — t z s A ) 2 2 
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Finally, we make use of the following Bessel function integral identity |)15f| : 

roo g 

/ ds — — — -J v (as)J v (bs) = I u (bc)K v (ac) (a > b), 
Jo s A + cr 

where I u and K v are modified Bessel functions. Taking, A 2 e 2 = 1 — A 2 , we have 

G{2N,2M,X) = ' 2 ^r D,2 {r') l - D,2 Ii^ 1 {^r l )Ki^ 1 ^r). (5.7) 



Apart from a multiplicative normalization constant, the expression in Eq. (|5.7|) is the 
final result for the Euclidean propagator. Let Q{r — > r') represent the spherically averaged 
amplitude for a free scalar particle of mass fi to propagate from some point on a sphere of 
radius r to some point on a sphere of radius r'. Note that this probability amplitude is not 
symmetric under the interchange of r and r'; when D > 1 it is more likely for a particle to 
propagate from a sphere of smaller radius to a sphere of larger radius than for the reverse 
to occur. This is because the final state of the particle propagating to the larger sphere has 
a higher entropy. This asymmetry does not occur in translationally invariant theories. Our 
final, properly normalized, result for the propagator is 

Q(r _> r>) = (r') D - 1 (rr') 1 - D / 2 In fl (fir < )KD f 2(fir > ), (5.8) 

where 

r > = max{r, r'} and r< = min{r, r'}. 



The normalization of the Green's function in Eq. (|5.8j ) will be verified in Sec. |VB . 

The propagation asymmetry in Eq. (|5.8|) is a continuum manifestation of the directional 
bias that is present in spherically symmetric random walks. Note that C n ^ m , the probability 
of walking from m to n [see Eq. fl2.1| )l, is not a symmetric function of m and n. Rather, it is 
the function v n _\ in Eq. ( |3.7| ) multiplying a symmetric function of m and n. The function 
v n -i represents the random- walk entropy associated with the volume of hyperspherical region 
n. The asymmetry in Eq. (p.8|) is a direct consequence of the asymmetry in C n ^ m . 
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B. Normalization of the two-point Green's function 



The free propagator in momentum space for a D-dimensional translationally symmetric 
scalar field theory is 



k 2 + n 2 

To obtain the coordinate-space propagator Q(r — r') we take the D-dimensional Fourier 
transform of the momentum-space propagator: 



2n) D J k 2 + fi 2 

W/V|r-r1)^i^(/i|r-r'|). (5.9) 



The coordinate-space propagator satisfies the Green's function differential equation 

(V 2 - /i 2 )£(r - r') = 5 (D) (r-r). 

Let us calculate the amplitude for a particle at r to propagate anywhere. We obtain this 
amplitude by integrating Q(r — r') with respect to r' over all space: 

f d D r'g(r-r') = - 2 . (5.10) 

We can now verify that Q{r — > r') in Eq. ( |5.8| ) is properly normalized by calculating the 
amplitude for a particle at radius r to propagate to any radius: 



poo 

[r — > r /A 



/•CO 

/ dr'Qii 
Jo 

= r l - D/2 KD^(nr) I dr' {r') D ' 2 In^(^r') + r 1 "^ 2 /^ (/ir) I" dr' (r') D / 2 Kn^(fj,r') 



r 

i 



7d (r)Ko-2 (r) + J d- 2 (rjifc (r) 



where we have the used the Wronskian identity for modified Bessel functions. This result 
agrees with that in Eq. ( |5.10|) . 
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C. Continuum derivation of spherically symmetric propagator 



In this subsection we derive the spherically symmetric propagator in Eq. (|5.8|) from the 
translationally symmetric propagator in Eq. ( |5.9|) by taking an angular average. To obtain 
the angular average we let 



|r — r'| = yr 2 + (r') 2 — 2rr' cos6. 
We then expand the modified Bessel function in Eq. (|5.9|) as a series in terms of Gegenbauer 



polynomials: 



' n=0 

If we then integrate over the angle 8, only the n = term in the series survives and we 
obtain the result in Eq. ( |5.8|) . The fact that we obtain the same two-point Green's function 
directly from our random walk model supports the validity of the uniform approximation 
for the probabilities P ou t(n) and P in {ji) in Eq. ( |3.2|) . 
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FIGURES 

FIG. 1. Comparison between P\ n (n) in Eq. (1.2) and the uniform approximation to Pi n (n) 
in Eq. (3.2) for D = 3 and D = 5. Note that the uniform approximation is exact at D = 1 and 
D = 2. 
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TABLES 



TABLE I. Actual and Predicted Values of J2n=i n Cn,t;m for t = 1000 and m = 1 



k D 



Actual 



Leading Behavior with 1st Correction with 2nd Correction 



1 
2 
3 
4 
5 



2.622502 x 10 1 

4.013810 x 10 1 

5.047526 x 10 1 

5.987220 x 10 1 

6.631717 x 10 1 



2.523133 x 10 1 

3.963327 x 10 1 

5.046265 x 10 1 

5.944991 x 10 1 

6.728353 x 10 1 



2.623133 x 10 1 

4.013327 x 10 1 

5.046265 x 10 1 

5.894991 x 10 1 

6.628353 x 10 1 



2.623133 x 10 1 

4.013823 x 10 1 

5.047527 x 10 1 

5.897220 x 10 1 

6.631718 x 10 1 



1 

2 
3 

4 
5 



1.051450 x 10 3 

2.040138 x 10 3 

3.001000 x 10 3 

3.943028 x 10 3 

4.870366 x 10 3 



1.000000 x 10 3 

2.000000 x 10 3 

3.000000 x 10 3 

4.000000 x 10 3 

5.000000 x 10 3 



1.050463 x 10 3 

2.039633 x 10 3 

3.000000 x 10 3 

3.940550 x 10 3 

4.865433 x 10 3 



1.051463 x 10 3 

2.040133 x 10 3 

3.001000 x 10 3 

3.943050 x 10 3 

4.870433 x 10 3 



1 5.352671 x 10 4 

2 1.219251 x 10 5 

3 2.019011 x 10 5 

4 2.914616 x 10 5 

5 3.892603 x 10 5 



5.046265 x 10 4 

1.188998 x 10 5 

2.018506 x 10 5 

2.972495 x 10 5 

4.037012 x 10 5 



5.346265 x 10 4 

1.218998 x 10 5 

2.018506 x 10 5 

2.912495 x 10 5 

3.887012 x 10 5 



5.352573 x 10 4 

1.219246 x 10 5 

2.019011 x 10 5 

2.914651 x 10 5 

3.892731 x 10 5 



1 3.205902 x 10 6 

2 8.237810 x 10 6 

3 1.500000 x 10 7 

4 2.342114 x 10 7 

5 3.344349 x 10 7 



3.000000 x 10 6 

8.000000 x 10 6 

1.500000 x 10 7 

2.400000 x 10 7 

3.500000 x 10 7 



3.201851 x 10 6 

8.237800 x 10 6 

1.500000 x 10 7 

2.340550 x 10 7 

3.338520 x 10 7 



3.205851 x 10 6 

8.237800 x 10 6 

1.500000 x 10 7 

2.342150 x 10 7 

3.344520 x 10 7 
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TABLE II. Actual and Predicted Values of ££°=i n k C nAm for t = 1000 and D = 2 





m 


Actual 




Leading Behavior 


with 1st Correction 


with 2nd Correction 


1 


1 


4.013810 


X 


10 1 


3.963327 


X 


10 1 


4.013327 


X 


10 1 


4.013823 


X 


10 1 




3 


4.019752 


X 


10 1 


3.963327 


X 


10 1 


4.013327 


X 


10 1 


4.019768 


X 


10 1 




5 


4.033597 


X 


10 1 


3.963327 


X 


10 1 


4.013327 


X 


10 1 


4.033639 


X 


10 1 




7 


4.055306 


X 


10 1 


3.963327 


X 


10 1 


4.013327 


X 


10 1 


4.055438 


X 


10 1 




9 


4.084814 


X 


10 1 


3.963327 


X 


10 1 


4.013327 


X 


10 1 


4.085163 


X 


10 1 


2 




2.040138 


X 


If) 3 


2.000000 


X 


10 3 


2.039633 


X 


10 3 


2.040133 


X 


10 3 




3 


2.046198 


X 


10 3 


2.000000 


X 


10 3 


2.039633 


X 


10 3 


2.046133 


X 


10 3 




5 


2.060336 


X 


10 3 


2.000000 


X 


10 3 


2.039633 


X 


10 3 


2.060133 


X 


10 3 




7 


2.082553 


X 


10 3 


2.000000 


X 


10 3 


2.039633 


X 


10 3 


2.082133 


X 


10 3 




9 


2.112848 


X 


10 3 


2.000000 


X 


10 3 


2.039633 


X 


10 3 


2.112133 


X 


10 3 
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